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Abstract 

We present a nonperturbative QCD calculation of elastic J ftp meson production in photon-proton 
scattering at high energies. Using light cone wave functions of the photon and vector mesons, 
and the framework of the model of the stochastic QCD vacuum, we calculate the differential and 
integrated elastic cross sections for 7 p — > J/ip p. With an energy dependence following the 
two-pomeron model we are able to give a consistent description of the integrated cross sections 
and the differential cross sections at low \t\ in the range from 20 GeV up to the highest HERA 
energies. We discuss different approaches to introduce saturation and find no specific effects up to 
energies presently available. We also calculate and compare to experiments the cross section for T 
photoproduction. 

PACS numbers: 12.38.Lg,13.60.Le 
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I. INTRODUCTION 

Photoproduction of J/if> mesons is at the borderline of soft and hard physics. On one 
side the mass of the charmed quark provides a relatively large scale, but on the other hand 
the size of the J /if) meson is determined not only by the mass of the charmed quark but also 
by the confinement mechanism which therefore cannot be neglected. Indeed if confinement 
effects could be totally neglected the size of the J / if) would be of the order of the Coulomb 
radius l/(m c a s ) which is considerably larger than the Compton-wave length l/m c . 

J /if) photoproduction has been treated extensively with methods of perturbative QCD 
0, |3], f|, H . In this paper we present a nonperturbative approach. This has the disadvantage of 
stronger model dependence, but the advantage that the process can be viewed from an unified 
point of view together with other processes already studied, and no use of external quantities 
like parton distributions is needed. The only intervening quantities are inherently calculated 
in the nonperturbative approach and no new free parameters have to be introduced. By 
comparing the successes and limitations of the perturbative and nonperturbative approaches, 
important insight in the transition region between the two QCD regimes can be obtained. 

The approach presented here is based on a functional integral treatment of high energy 
scattering || where the functional integrals are evaluated in an extension of the stochastic 
vacuum model |7|, §. The method has been applied with great success to calculate differential 
and total cross sections for many processes. 

Although the size of the T is much smaller than that of the J /if) and hard contributions 
are expected to be dominant, we nevertheless also calculate photoproduction of the T-meson 
in our model, obtaining reasonable agreement with experiment. 

Our paper is organized as follows. In Sect. 2 we discuss shortly the main features of the 
underlying nonperturbative model and present the final formulae for dipole-dipole scattering, 
which is the basic ingredient in our approach. We also give convenient parametrisations of 
our theoretical results and discuss the inherent limitations of the model. Our numerical 
results and a comparison with experiment are given in Sect. 3. The paper closes with the 
discussion in Sect. 4. In the Appendix we collect some useful formulae concerning wave 
functions. 

II. THE MODEL 

A. Nonperturbative treatment of scattering amplitudes 

In this subsection we give a short review of the main ideas behind our nonperturbative 
model for soft high energy reactions and present the final formulae. For more details we refer 
to the original literature || [| and reviews [JR], [TT|, [RJ] . The functional integral approach to 
soft high energy scattering |J starts from the scattering of a highly energetic quark in an 
external colour field B M (x). Along its path T the quark picks up the non-Abelian phase 

—in f Bcix 
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where the expression is path ordered. Here and in the following we express through bold 
face letters matrix valued quantities like 

B M (s) = £ \\ C B°{x) , (2.1) 
c=i A 

where Ac represents the Gell-Mann matrices. 

According to the functional integral approach to quantisation the scattering amplitude 
of two quarks can be obtained by averaging these phase factors of two quarks with the 
exponential of the action as weight. Formally this can written as the functional integral 

J VBe~ i9 ^ Bd V^/r 2 Bdx exp[-iS QCD } = (e~ i9 ^ B< V* 9 4 B ^ )fi . (2 . 2 ) 

This quark-quark-scattering amplitude is neither observable nor gauge invariant. In order 
to construct a gauge invariant expression we consider the scattering of two colour neutral 
quark- ant iquark states, so called colour dipoles. This leads to the expectation value of two 
Wegner- Wilson loops ||, as depicted in Fig. 




FIG. 1: The scattering of two dipoles. 



The four corners of the loops have the coordinates 

d : (-T,-T,x 1 ),(T,T,x 1 ),(-T,-T,x' 1 ),(T,T,x' 1 ) (2.3) 
C 2 : (-T,T,x 2 ),(T,-T,x 2 ),(-T,T,x' 2 ),(T,-T,x' 2 ) (2.4) 

where the arrow always indicates vectors in transverse space and T goes to infinity. 
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The relative and centre coordinates are introduced as 



Ri x i x $ , X{ x i ~\~ z% R% . (2.5) 

— * 

The vector Ri denotes the transverse extension of the loop i and the quantity Zj with 
< Zi < 1 will later be identified with the longitudinal momentum fraction of the quark. 
The impact parameter vector h is defined by 

b = X 1 -X 2 . (2.6) 

With this definition the t-dependent scattering amplitude can be obtained as the two- 
dimensional Fourier transform with respect to the impact parameter |L3] . 



The basic element of the scattering matrix for colour- singlet quark-antiquark dipoles is 
the expectation value of two loops, 

S (b, R,, z u R 2 , z 2 )) - i {w[Ci] y {w[C2]) , (2-7) 



with 

W\Ci] = tr Pexp\-ig 

is, 



W[Ci) = tr P exp[-ig I Bdx] (2.8) 



where Si the border line of loop i. 

We will discuss later how we pass from the dipole-dipole amplitudes to hadronic scattering 
(or photoproduction) amplitudes by integrating over light-cone wave functions. 

The expectation value of the loops is approximately calculated using an extension of 
the stochastic vacuum model (SVM). In this model it is assumed that the long-distance 
behaviour of QCD can be approximated by a Gaussian stochastic process with the gluon 
field strength as the stochastic variable. This model yields confinement in non-Abelian 
gauge theories and is in conformity with the Mandelstam-t'Hooft 15 picture of string 
formation through monopole condensation; we refer to flJJ] for a detailed review. In order to 



apply the model for the evaluation of two loops it has to be extended, since in a stochastic 
process with non-commuting variables the higher cumulants which must vanish in a Gaussian 
process are not uniquely defined. 

In order to pass from the gluon potential occuring in the line integrals in Eq. ( [2.8|) to 
the gluon field strength we have to apply the non-Abelian Stokes theorem. This implies a 
special choice of a surface containing both loops as borders. This choice is not unique, and 



in this paper we use one proposed in 0]; for other possible choices, see [0. 

The general form of the basic correlator which determines the full Gaussian process is in 
our approach given by M 



(:^W(e" i9/: Adz ) C i£0O:> 

d A k 

W) 



8 ab 8 cd (g 2 FF) B J ^- e-^~ x '\{g» P g v „ - g^g vp ) « iD{k 2 



'9v<jk^kp + gupk^ku — g^pkyk^ + g^kykp^iX — k) i- 



dk 2 



(2.9) 
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For the correlation functions we take the form proposed in || 



D(k 2 ) 



277T 4 fc 2 

9tt 4 

2« 2 (^ 2 + i&) 3 



with the parameters 



0.346 fm 



(g 2 FF)a i = 23.5 



0.74 



(2.10) 



(2.11) 



where a is the correlation length, which are in agreement with lattice results [|T^, [H]]. This 
same parameter set has been used for many applications of the model to hadron-hadron, 
photon-hadron and photon-photon high energy interactions and will be used throughout this 
paper. 

The formalism set up above allows us to calculate the scattering matrix in terms of the 
correlator (2.9) using the assumptions of the stochastic- vacuum model. The most straight- 
forward way is to expand the exponentials occuring in the Wegner- Wilson loops of Eq. Q2.7| ) . 
We then obtain 



5(6,1,2 
x 
+ ... 



| i tr(A a A fe )tr(A c A^ 
do KX F c ; 



144 



Si 



s 2 



Si 



da^Fl 



da KX F: 



s 2 



B 



(2.12) 



where the dots represent products of more than four field-strength tensors and 1,2 stands 
for Ri, zi, i?2, z%. Inserting the expressions ( |2.10|) and performing the surface integrals we 
are finally lead to 

S(M,2) = 1-| X 2 (M,2), 

with 



(2.13) 



X (M,2) = -(g 2 FF)(l(x 1 ,x 2 ) + I(x' 1 ,x' 2 

- /(fi,4) - 

Using the special form of the correlators in Eq. (|2.10|) we obtain 



(2.14) 



I(xi,x 2 ) = |vtk J dv (\vx 2 — Xi\ 2 K 2 (\ 1 \vr 2 — xi\) 

+ \x 2 - vx 1 \ 2 K 2 (X~ \X2 - VXi\ 
+ (1 — K)7rA 2 |a;2 — xi| 2 i^3(A _1 |x2 — xi|) 



(2.15) 



where A = (37r/8)a , and K 2 and K 3 are modified Bessel functions. 

A more refined method to treat the two traces has been developed by Berger and Nacht- 
pOj . The main idea is to interpret the product of the two separate traces in Eq. (|2.12|) 



mann 



over 3x3 matrices (generators of SU (3) in the fundamental representation) as one trace 
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Tr 2 in the product space of the two fundamental representations of 577(3); that is Tr 2 acts 
in SU(3)®SU(3). Thus 

5(6,1,2) = |Tr 2 (exp(-^/ 5i rfa^(|A a ®l)) 

x exp ( - ig i da^F^l ® iA c ))) fi . (2.16) 
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The two exponentials commute in the product space and we can write the right-hand side 
of this equation as a single exponential 

5(6,1,2) = |Tr 2 (exp(-z(?/ 5 rfa^F^)) B , (2.17) 

where the surface integral extends over the surfaces Si and S 2 , and takes its values in 
the product algebra of SU(3)<S>SU(3). We can now make a cluster expansion with stochastic 
variables from the product algebra and using the correlator ( [2 .ICQ we finally obtains 



5(6,l,2) = fe-3 i * + ie§ ix , (2.18) 



where \ = x(6, 1,2) is the same as given in Eq. (|2.14 ). 



In the following we refer to the first method leading to Eq. ( |2.13| ) as the expansion 
method, and the second method leading to Eq. ( p. 18 ) as the matrix-cumulant method. 



In a very loose sense we can consider the quantity x as representing the exchange of 
a nonperturbative gluon. In this sense the expansion method takes only into account the 
exchange of two nonperturbative gluons. The matrix cumulant method then takes also 
into account multiple gluon exchange; it will automatically satisfy unitarity constraints for 
hadronic cross sections. An expansion of Eq. (|2.18| ) in \ yields in leading order the result 
obtained in Eq. ( |2.13| ) above, namely 

5(6,1,2) = l-| x 2 -... (2.19) 

We also consider the eikonal unitarisation method 

5(6,1,2) = e~ x2 / 9 , (2.20) 

which is similar to forms used to study saturation effects 0, ^T[]. Like the matrix-cumulant 
method it corresponds to multiple exchange, but here the exchanged objects represented 
by x are coupled in such a way that always a pair forms a color singlet whereas in the 
matrix-cumulant method only the entirety of the exchanged objects has to form a colour 
singlet. The three different methods are illustrated in Fig. |2| . 

The matrix element for dipole-dipole scattering with momentum transfer q is then given 

by 

T fi (s,t,R h zi,R 2 ,z 2 ) = -2%s J d 2 be iq ^(S(b, 1,2) - l) , (2.21) 

with t = —q 2 . 

The treatment of three quarks in a colour- singlet state is analogous to the scattering 
of two dipoles, but technically more involved. We refer to the literature || |2(| for the 
corresponding results. 

If the parameters of the stochastic-vacuum model are taken to be independent of the 
energy, the resulting cross sections turn out to be independent of the scattering energy too. 
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a) 



b) 



c) 



FIG. 2: Illustration of a) the expansion metod, b) the matrix-cumulant method and c) the eikonal 
method. The gluon line with a x corresponds to the expression ( [2.14 ) and can be viewed loosely 
as the exchange of a nonperturbative gluon. 



The observed energy dependence has therefore to be introduced by hand. One way is to 
make the size of the hadrons energy dependent || ^2[. This has been shown to be very 
simple and effective for purely hadronic processes. We here adopt the two-Pomeron approach 



of Donnachie and Landshoff [ 25fl , coupling a soft Pomeron with intercept 1.08 to large and 
a hard Pomeron with intercept 1.42 to small dipoles G41 125 . 



Specifically we introduce an energy dependence of the quantity Xi Eq. Q2.14-D 

/ \ W 2 



x(M,2) 
x(M,2) 
x(M,2) 



X h (b, l,2,s) = X (b, 1, 2) ( - 1 for R x and R 2 > r c 



X E (6,1,2, S ) = X (6,1,2) 



sR[ 



2\£h/2 



X E (b,l^s) = X (b,l,2)( sRlR2 



\ S r c 



for Ri < r c , i = 1 or 2 
| for Ri and R 2 < r c 



(2.22) 



with e s = 0.08, = 0.42 taken from 23] and r c =0.22 fm taken from a treatment of the 



proton structure function [25|. 



B. Wave functions and hadronic reactions 



We are interested in reaction amplitudes where the external particles are physical hadrons 
and photons. We obtain such amplitudes from the dipole-dipole scattering amplitude Q2.21 



by integrating over all dipoles sizes with the light-cone wave functions of the participating 
particles as weights. 

A meson or a photon is here described by a light-cone wave function i[) n (Ri, Zj) of a quark 
and an antiquark with relative transverse coordinates Ri and quark longitudinal momentum 
fraction Zi. The meson or photon scattering amplitude for the reaction a b — > c d is then 



obtained from the dipole-dipole scattering amplitude T/i(s, t, Ri, z±, R 2 , z 2 ) of Eq. ( [2.2 1|) by 



T ab ^ cd (s, t) = -2is J d 2 be^ rb (S ab ^ cd (b, s)-l), (2.23) 



where 



S ab ^ cd (b,s) = J d 2 R x J d 2 R 2 dz i J Q dz 2 VCC^i, zi)i>a{Ri, z\) 
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(2.24) 



The normalization is such that 

|pT^ |r " (M)|2 - ^ 

For the heavy quark content of a photon the perturbative expressions for the wave func- 
tions are reliable since the charm quark mass sets the scale. Although we need in the present 
work only with the transverse wave function, we give for completeness the expressions for 
both transverse and longitudinal photons, the latter being needed for electroproduction of 
mesons. For photons of helicities 1, -1 and 0, we write respectively 



2tt 



iee 1 (z8 hj+ 8 h _ - z8 h _8- h ^)Ki{er) + m f 8 hy +8- h>+ K G (er) 



(2.26) 



if) 7 -i(Q 2 ;z,r,9) 



2tt 



iee t6 (z5 hi+ 5- h ^ - z8 h ^8j l ^)K x {er) + mf8 h -8j l _KQ{er) 



and 



where 



'3a 



2tt 



Izz) S h ,-h Q ^o(er) 



zzQ 2 + m 2 



f ' 



^2.27) 



(2.28) 



^2.29) 



and mj is the quark mass and if is the quark charge in units of the elementary charge for 
each flavour /; K , K\ are the modified Bessel functions. 

Meson wave functions are more model dependent than photon wave functions and espe- 
cially the spin structure can be quite complicated p6| . In this paper we take for the vector 
mesons the spin structure from the vector current leading to similar expressions as for the 
photon JIB], namely 



z,r) 



ipv,-i{z,r 



VU, r)\ iu 2 re l9 (z5 h 



-5h - z5 h -5 



) v {z,r)\iuj 2 re ^(zS^+Si 



h 

z8 h ,-8 



+ m f 8 h ,+8 hj+ ) 
.) +mf8 h -8j l 



and 



tpv,a(z,r) = (f)(z,r) v {ujAzz8, 



h,-h 



[2.30) 



(2.31) 



Here ±1 and denote transverse and longitudinal polarizations of the vector meson, and h 
and h represent the helicities of quark and antiquark respectively. 

The functions <py(z,r) are constrained by the normalisation condition and by the elec- 
tromagnetic decay width, as described in the Appendix. Making for the r dependence a 
Gaussian ansatz, the parameters are then completely determined by the two conditions just 



mentioned. For the z dependence we make two ansatze: one is suggested by the phenomeno- 
logically very successful! Bauer-Stech-Wirbel model [27 



N 



z(l — z) exp 



2lu 2< 



exp[— -uj 2 r 2 ] 



(2.32) 



The other choice is obtained in the spirit of Brodsky-Lepage |28| from a non-relativistic wave 
function in the rest frame where the transition to the light-cone coordinates is achieved by 
the replacement of the relative momentum k of the two constituents by the transverse 
momentum kx and the longitudinal momentum fraction of the quark according to 



k 



k% + m 



1 



4z(l - z) 



m 



f 



(2.33) 



This procedure only makes sense for a finite quark mass and leads to a more photon-like z 
dependence of the wave function. We then write 



N 



>BLU,r) 



exp 



m 2 Az 



1)2 

2/ 



2uo 2 z(l 



exp[— 2z{l — z)uj r ] 



(2.34) 



My and rrif represent respectively the vector meson mass and the the quark mass. 

For the proton we should use a three quark wave function which in principle poses no 
problems |5J. Earlier investigations have shown however that a quark diquark structure of 
the proton leads to phenomenologically consistent results and most applications have been 
made in this picture. The quark-diquark nucleon wave function can be treated like a quark- 
antiquark meson wave function and again we make a Gaussian ansatz for the wave function 
and fix the longitudinal momentum fraction at z ~ 1/2 , namely 



MR) 



]_}_ e -r*/{2S p ? 

2tt St, 



The transverse size parameter S p of the proton is choosen to be 

S„ = 0.74 fm , 



(2.35) 



(2.36) 



which leads to a good description of pp scattering and also to a good proton form factor 
I- 

For the charm mass we adopt the value used in a previous analysis of structure func- 

(2.37) 



tions p5| , 



1.25 GeV . 



This value lies in the middle of the range of masses of the modified minimal subtraction 
scheme PD[. The dependence of the charm mass is discussed in the Appendix. 

This finishes the description of the general model, with characterization of all parameters. 
It should be noted that in the original application of the matrix cumulant method a 
slightly different set of parameters was used in order to get an optimal overall description of 
the elastic proton-proton cross section. Showing in the present paper the values of amplitudes 
and cross sections evaluated with only one set of parameters we wish to emphasize the 
magnitude of the influence of multiple scattering mechanisms. 
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C. Photoproduction of vector mesons 



Using the results of the two previous sub-sections we can calculate the scattering ampli- 
tudes for photoproduction of vector mesons. 

Inserting the wave functions of Eqs . (|2 . 2 7| ) , ( |2 . 30] ) and (|2.35|) into the production amplitude 
( |2.24| ) we obtain 

S^ Vp (b) = J d 2 R 1 J dz x J d 2 R 2 p^ x (z 1 ,R 1 )\i; p (R 2 )\ 2 S(b,z 1 ,R 1 , 1/2, R 2 ) , (2.38) 
where p 7 y t \(zx, R^) is the photon-vector-meson overlap function 



P-f,v,x(zi, R±) — e q — — 4>bsw/bl{z-,Ri 

Z7T 



x [e u] l R x \z l + (1 - zy\K 1 {m s R x ) + mjK (m f R^j . (2.39) 

The overlap densities are independent of the angles 9i and 9 2 of R\ and R 2 , and the 
functions x(6, 1, 2) change sign if R\ or R 2 is reversed, as can be seen easily from Eq. (|2.14|) 
by noting that Ri — ► — Ri corresponds to (xj,x^) — * (x'^Xi). Inserting the result (|2.18|) of 
the matrix cumulant method into Eq. ( |2.38| ) we see that after integration over the angles 9\ 
and 9 2 only the even terms in x survive and the exponentials in Eq. ( |2.18j ) can be replaced 
by cosines. Therefore we may insert in Eq. ( |2.38j ), as result of the matrix cumulant method, 

5(6, 1,2) = | cos + | cos (fx) = 1 + | [ cos + 2] [ cos (| X ) - l] (2.40) 

This form guarantees that the amplitude remains inside the unitarity bounds if x becomes 
large. Differences with respect to the expansion method, Eq. ( |2.13|) , increase at high 
energies. 

If the energy dependent expressions in Eq. ( |2.22|) are introduced it turns out that the 
hard part of the proton, R± < r c gives for energies below the TeV region only a negligible 
contribution, so we have only to consider the two cases R\ < r c and R\ > r c . 

Now all formulae are set up and all parameters are fixed. 

Our results for the differential and total cross section of J/ip photoproduction in the 
expansion method are written 

wr (Mt)(-) +*(«)(-) ) . (2.41) 

The soft and hard amplitudes for the case of the BSW wave function can be conveniently 
expressed by the parametrisations 

10.82 exp(-1.07|t|) 



A s (t) 



1 + 3.391*1 



Mt) = 2 - ^y) , (BSW) , (2.42) 

where t is in GeV 2 and A St h{t) in v^nb/GeV, sq = (20 GeV) 2 . For the case of the BL wave 
function the forms are very similar 

11.17 exp(-1.06|t|) 



AM = 

w l + 3.26|t| 
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exp. 


matrix-cum. 


■1 l 

eikonal 


w 


A 


A 


A 


GeV 


nb/GeV 2 


nb/GeV 2 


nb/GeV 2 


20 


167.6 


153.7 


141.1 


200 


896.4 


790.2 


693.3 


1000 


5657 


3983 


3008 



TABLE I: Result for A in the parametrisation ( 2.44 ) for the expansion, the matrix-cumulant and 
the eikonal methods. 



These very convenient parametrisations reproduce the exact results with an accuracy 
always better than 5 percent. We see that the t-dependence of the differential cross section 
predicted by our model is not a pure exponential, showing a curvature in a logarithimic 
scale. Since the parameters expressing the t-dependences of the hard and soft parts of the 
amplitude have rather similar values, the shape of the angular distribution (and the slope 
parameter) depend only weakly on the energy (we recall that we work with small values of 
\t\). The values of A s (t) and Ah(t) are shown in Fig. ([3D, in solid line for the BWS and 
in dashed line for the BL wave function. These forms, together with Eq. (|2.41|) contain all 
results for the expansion method that will be used for comparison with experiments. At low 
energies the soft contribution is several times stronger than the hard one. The hard part 
reaches the soft part for W about 240 GeV. 

Since the differences between the two kinds of J/ip wave functions are very small, from 
now on in the present paper, we will use only one of them, namely the BSW wave function. 

The results for the matrix cumulant method [20] cannot be parametrized so easily, due to 
the nonlinear dependence of the amplitudes on the quantity x E ip, 1> 2, s), Eq. ( |2.22| ). That 
is, we cannot factorize the s and t dependences of the soft and hard parts, as in Eq. fl2.41| ) 



The same is true of the amplitudes obtained with the eikonal form of Eq. (|2.2U| ). In Fig. 
|3] we draw also the amplitudes for these two unitarization procedures, for a fixed energy 
W = 20 GeV, where we see that the form of the t-dependence (in a limited \t\ interval) is 
not much affected by the saturation corrections. The same is true for higher energies. 
Accurate parametrisations of the differential cross sections can be written in the form 

da_ _ Aexp(-a\t\) . , 

d\t\~ l + b\t\ ' 1 ] 

yielding results summarised in Tables |, ||. We found that parametric forms with dipole 
factor like 1/(1 + b\t\) 2 do not lead to equally accurate representations for da/d\t\. 
For the integrated production cross section we obtain with the expansion procedure 

/ S \0.84 / S \0.16 / S \0.50 

ct = 0.552( — ) + 13.539( — ) +5.465( — ) . 2.45 
V V V V v s y 

In the upper part of Fig. f| we show the results for the integrated cross sections using the 
same input-parameter set of Eqs. ( |2.11 ) and ( |2.36| ) for the expansion method of Eq. fl2.13p 



and for the matrix cumulant expression of Eq. ( |2.40| ) and for the eikonal form of Eq. Q2.20D . 
The curves show the influence of possible multiple scattering corrections. In the lower part 
of Fig. [| we instead choose slightly different values for the gluon condensate: (g 2 FF)a 4 = 
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expansion 


matrix-cumulant 


eikonal 


w 


a 


b 


a 


b 


a 


b 


GeV 


GeV" 2 


GeV" 2 


GeV" 2 


GeV- 2 


GeV" 2 


GeV" 2 


20 


2.72 


9.0 


2.83 


10 


3.00 


10 


200 


2.62 


9.2 


2.87 


10 


3.15 


10 


1000 


2.56 


9.3 


3.38 


10 


3.78 


10 




FIG. 3: i-dependence of soft and hard production amplitudes at the energy W = \fs = 20 GeV. 
The lines can be parametrised in forms A(0) exp(— a|i|)/(l + b\t\), as explained in the text. The 
solid and the dashed lines are the results for the expansion method Eq.( 2.13 ), using respectively 
BSW and BL wave functions for J /if). The dotted line corresponds to the matrix-cumulant method 
Eq.(2.40) and the dot-dashed line to the eikonal method Eq.(2.20), both using BSW wave function. 



25.0 and 26.0 respectively for Eqs. ( |2.40| ) and ( |2.20| ). This last figure shows that in the 
present experimental range, below 300 GeV, unitarity corrections are easily compensated by 
slight changes of one single parameter. Small changes in the value of the charm quark mass 
that enters in the wave function has important effects of the same kind. 

We may conclude that these saturation effects are not clearly observed in elastic photo- 
production of J /if) for energies up to 300 GeV. 

It is interesting to investigate the behaviour of the integrated cross section at higher 
energies. While with the expansion method the cross section increase according to Eq. 
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FIG. 4: Integrated elastic production 7 p — > J /tp p cross section. In full line the result with the 
expansion method, Eq. fl2.13 ). Dotted and dash-dotted lines correspond respectively to the matrix 
cumulant method of Eq. ( |2.40| ) and to the exponentiation form of Eq. fl2.2tip . In the upper figure 
the same (g 2 FF)a = 23.5 for the gluon condensate is used for the three curves. In the lower 
figure we show that the unitarization effects are compensated, in the present experimental range 
< 300 GeV, by slight changes of parameters: {g 2 FF)a i = 25.0 and 26.0 respectively for Eqs. 
(|2~40| ) and (^20|) . The behaviour at high energies is explained in the text. 



( p.45| ), in the unitarized cases the behaviour in the whole range 20 < W < 1000 GeV is 
described by the slower behaviour 

a(W) = a(W ) x [1 + C(^-) 5 log , (2.46) 

with 5 = 0.7 for the matrix cumulant (dotted lines) cases and = 0.6 for the exponentiation 
procedure (dot-dashed lines). The values of the multiplicative constant C vary a little with 
the choice of {g 2 FF)a i = as follows: C=0.36 and 0.35 in the matrix cumulant method respec- 
tively for 23.5 and 25.0 , and C=0.43 and 0.41 in the exponentiation procedure respectively 
for gluon condensates 23.5 and 26.0. 

If we restrict ourselves to the experimental energy range from 20 to 300 GeV, the results 
obtained with the expansion method can be represented with an effective power like form 

<t(W) = 14.6 + 9.37(— ) L08 . (2.47) 



The functional integrals occuring in the basic expression (|2.71) for the S-matrix are in our 
approach approximately evaluated with the help of the stochastic vacuum model. This model 
is an approximation appropriate only for the soft part of QCD. The higher the momentum 
transfer the more important will become contributions of hard gluons. Therefore we expect 
our calculations to be best at small momentum transfer. Of course we cannot predict 
the exact scale where hard scattering becomes important but investigations of the gluon 
distributions in hadrons and virtual photons |31] indicate that the hard component becomes 



as important as the soft part at a transverse momentum of the gluons \kx\ ~ 1 GeV, so that 

a safe limit for soft process is y\t\ <C 1 GeV. We also expect that for photoproduction of 
T-mesons even at small momentum transfer hard gluons play a more important role . There 
the small size of the T suppresses the contribution of the soft (nonperturbative) gluons as 
compared to the hard ones. 

Another limitation of our approach is given by the WKB approximation underlying the 
nonperturbative approach to scattering. Here the path is assumed to be the classical 
one, being nearly a straight line. This also determines the momentum transfer below which 
the model can be safely applied. The WKB approximation should be more reliable at high 
energies, and this restricts the application of the model to energies higher than -y/i « 20 
GeV. 

Apart from the approximations characteristic of the stochastic vacuum model, also the 
wave functions are treated in a very simplified form. From previous experience we expect 
the reliabillity of the model to be about 10 % in the amplitudes. Thus, putting all this 
together, discrepancies between theory and experiment below 20 % for the production cross 
sections could be considered as natural. 

III. COMPARISON WITH EXPERIMENT 

The experimental |32] forward production cross section (da/dt)(t = 0) is shown in Fig. 



II Only the HI Collaboration measures directly this quantity, represented by their < 
\t\ < 0.053 GeV 2 bin . The values reported by Zeus [|33| correspond to the extrapolations to 



\t\ = of their fitted straight lines, as will be discussed later (Fig. 
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FIG. 5: Forward differential cross section (da/dt)(t = 0) for the reaction 7 p — > J/ip p. The 
line represents our results with the expansion method. The experimental results are from the HI 



Collaboration [32| for the momentum transfer in the interval < \t\ = 0.053 GeV 



A direct comparison with both ZEUS and HI data in the forward direction can be made 
at \t\ = 0.1 GeV 2 , and this is presented in Fig. || Here we see a very satisfactory agreement 
between theory and ZEUS and HI data, with a remarkable theoretical description of the 
energy dependence in the whole range from 40 to 260 GeV. 

The integrated production cross section is presented in Fig. ^|, showing again a reasonable 
agreement of our model with all data up to nearly 300 GeV. The solid line corresponds to 
our calculation using the \ 2 expansion of Eq. ( |2.13|) . As we have already shown in Fig. |], 
the effects of saturation introduced by the matrix cumulant or exponentiation method are 
not very large, are within our expected errors, and can be compensated by slight change in 
parameters. In the figure the dashed line represents the results obtained with the matrix 
cumulant, using the standard set of parameters. We have shown in Fig. |] that choosing 
(g 2 FF)a 4 = 25.0 , instead of 23.5, the result of the matrix-cumulant coincides in this energy 
range nearly with the solid line in Fig. |7|. To show the influence of the value of the charm 
mass in the calculation, this figure includes also (dot-dashed line) the results obtained using 
a value m c = 1.2 GeV in the wave functions. 

The HI and ZEUS data on integrated elastic cross sections put together can be fitted by 
either of the two forms (in nb) 



<r(W) = 25.8 



or 



W\o.7i 



W\ .43 - ,W 



a(W) = 26.0 + 17.14 (_)- log(-), 
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FIG. 6: Differential cross section (da / dt){\t\ = 0.1) for the elastic production 7 p — > J/-0 p 
at |t| =0.1 GeV 2 . The solid line represents our calculation with the expansion method. The 
experimental results are from ZEUS [33] and HI (3^] collaborations. Zeus data are separately 
reported for observations of J/ip — > and J/tp ^ e + e~ decays. 



with the same deviation \ 2 — 1-73 per degree of freedom. 

In Fig. |8| we compare our calculations of differential cross sections with the most recent 
data from HI and ZEUS experiments. The energy range is from 40 to 260 GeV, distributed 
in 9 bins, and the \t\ range is (0 — 1.6)GeV 2 . Where data from different experiments in 
similar energy bins are available, we plot them together in the same figure. We emphasize 
that our calculation contains no adjustable parameters and that these measurements cover 
a wide energy range. 

As a quite general feature we see that for small values of the momentum transfer t the 
agreement between experiment and theory is quite satifactory for all energies. For larger 
values of \t\ agreement with HI data is still satisfactory, but at the lower and middle energies 
our results are in general below the new ZEUS data. At the highest energies, the agreement 
improves. Since in our model we take into account only nonperturbative effects it is not 
surprising that at larger momentum transfers where harder gluons are exchanged, some 
contributions are missing in our calculation. In p|J a good fit up to \t\ « 6 GeV 2 is indeed 
obtained. 

Zeus gives exponential fits in t to their data, which are indicated by dotted lines in our 
plots in Fig. (P) . From this fit the Zeus collaboration obtains the forward differential 
cross section and then evaluates the integrated elastic cross section. The plots show remark- 
able differences at \t\ = between our results (with the epansion method) and the ZEUS 
extrapolated values obtained with a straight line. 

As we have shown in the general description of the calculation, the differences between 
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FIG. 7: Integrated cross section for the reaction 7 p — > J/ip p. The experimental data are from the 
HI |0] and Zeus collaborations. The point marked with a star at IF = 16.8 GeV represents 



results from fixed target experiments [35, [3J], 37]. The solid line is obtained with the expansion 
method. The dashed line shows our results using the matrix cumulant method, with the same set 
of parameters, as in the top part of Fig. |||. To show the influence of a different choice for the charm 
mass, the dot-dashed line gives the results obtained with the expansion method, using m c = 1.2 
GeV. 



the expansion method and the matrix cumulant and exponentiation methods are not very 
large and could be absorbed in a slightly different choice of parameters. We do not include 
the corresponding lines in Fig. [8] in order not to overload the plots. 

The data presented above are the most recent HERA data on J/ip photoproduction. 
Some of the pioneering fixed target experiments of 20 years ago were made at energies near 
t/s = 20 GeV , which is at the border of the range appropriate for our calculations. We 
have included in Fig. ^ a point at about 16.8 GeV representing this experimental effort 
35|, |36|, [37[|, showing that it fits well in the sequence of higher energy data points. As a 



historical tribute, we show in Fig. ^| the \t [-dependence of the differential cross sections 
obtained in some of these experiments. In this figure we draw together our curves for the 
expansion method of Eq. ( |2.13| ) (solid line), for the matrix cumulant result of Berger and 
Nachtmann [5D| given by Eq. fl2.40p and for the exponentiation procedure of Eq. ( |2.20| ). 

We have also applied the model to photoproduction of T (IS) mesons. As expected the 
results of our calculation are below the central values of the experimental data from Zeus 
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FIG. 8: i-dependence of differential cross sections for the reaction 7 p — > J/ijj p. Data are from HI 
p2| and ZEUS [33] experiments. The solid curves represent our calculations, as described in the 
text. The dashed straight lines are Zeus fits in^ach of their energy bins. 
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FIG. 9: i-dependence of differential cross sections for the reaction 7 p —* J ftp p at the low energies 

~~ I, |37l]. The solid, dashed and dotted lines represent our 



of some fixed target experiments [35 



results using Eqs. (EDI), (Egg) and Eq. |L20|) respectively. 



| 3q| and HI 32| , but the large errors do not allow a definite conclusion about the need 
of other contributions. We find it nevertheless quite astonishing that the nonperturbative 
model seems to yield at least a substantial fraction of the cross section for such a hard 



process. Our results shown in Fig. [10] can be compared with perturbative calculations |39 



40| of the same process. In [3l| the LO calculation taking into account the skewed parton 



distrubution in some way obtains results on the lower edge of the error bars, while in |4(^ 
additional contributions lead to values just through the central values of the measurements. 



IV. SUMMARY AND DISCUSSION 

The main purpose of this paper is to investigate to what extent photoproduction of 
J/ijj(lS) mesons can be described in a purely nonperturbative QCD model. The model 
has been tested before, and the parameters used in the present work have been taken from 
previous publications on different processes. They were determined mainly from results of 
QCD lattice calculations and hadronic properties. Our parameter free results are compared 
with the more recent experimental data obtained at HERA by the HI [32] and Zeus [33| 



collaborations. Comparison of the calculated values of the differential cross sections with 
the data in nearly forward directions is presented in Figs. [| and |6|. The values reported 
by Zeus at |t| = are not direct measurements, but rather extrapoled values from linear 
fits, as shown in Fig. ^. In our view, these linear extrapolations do not account for possible 
structure in the very forward direction, and may lead to underestimated values, so that Fig. 
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FIG. 10: Integrated cross section for the reaction 7 p — > T p. Solid : our result with mj = 4.2 
GeV, dashed with = 4.4 GeV. Experimental points from Zeus p3] and HI [32] 



[5] only shows the HI data. At \t\ = 0.1 GeV 2 both HI and Zeus direct measurements exist, 
and Fig. |^ shows that the theoretical model gives excellent description of the magnitude 
and of the energy dependence at low \t\. 

Fig. |S] shows that at larger \t\ our calculations agree well with the HI measurements and 
are below the Zeus points for energies up to W = 120 GeV. At the highest energies the 
agreement with Zeus data is good, except that our model predicts a curvature, with a rise 
in the very forward direction. A clarification of this point is very important, since many 
models can only calculate the forward scattering amplitude and the test of the corresponding 
theoretical ideas depend on this kind of experimental data. Our model leads to a peculiar 
form-factor dependence which is exhibited in Fig. ^ and represented by the parametrizations 
of Eqs. ( |2.41|) , ( |2.42| ) and ( |2.43| ). The t-dependence of the differential cross sections shows 
similar curvatures, represented by Eq. (|2.44|) . 

We recall that our model in meant to be tested in the low \t\ range, say below 1.0 GeV 2 . 
For larger momentum transfers the differential cross section values are already 100 times 
smaller than in the soft region, and an additional genuinely hard contribution, however small 
(without influence on the integrated cross section), may play here an important role. 

The integrated cross section is in reasonable agreement with the experimental data, as 
shown in Fig. [7]. 

We have studied particularly the influence of saturation. We have several methods based 
on the stochastic vacuum model to evaluate our model. 

• The expansion method which corresponds loosely speaking to an exchange of two 
,,nonperturbative" gluons. 
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The matrix cumulant method p(J which takes in some way multigluon exchange into 



account and respects the unitarity constraints for hadronic cross sections. 
The usual exponentiation of the profile function, as in Eq. ( [2.20 ). 



The methods are illustrated in Fig. ^] and compared in Figs. [3]and|4]. The energy dependence 



was introduced based on the two-pomeron model of Donnachie and Landshoff |2"3fl . In the 
expansion method it leads to a power like increase of the integrated cross section, while the 
matrix cumulant and in the exponentiation methods there is saturation due to the unitarity 
constraints on the dipole cross sections, inherently respected by the approach. Up to energies 
about 1000 GeV, the energy dependence in these unitarity controlled calculations is very 
well parametrized in forms given by Eq. ( |2.46| ). In Fig. |7| our calculations are compared 
with the data. 

It should be noted that in photoproduction processes the Froissart theorem cannot be 
proved and that a powerlike increase does not contradict fundamental theorems of local 
quantum field theory. 

At present energies and with present accuracies the data show no indication of saturation 
and are well compatible with the powerlike behaviour of the unconstrained two pomeron 
approach. But a comparison with the theoretical approaches that introduce saturation 
effects shows that these are also compatible with the experiment, considered some 20 % 
allowed variation in our model. Rather small changes of parameter values may account for 
differences among calculation procedures. 

We have also tested two different types of wave functions and found only small differences 
in relevant results. 

We conclude that, in the energy range available at HERA, the methods here studied yield 
results compatible with experiment. As can be seen from Fig. [|, sizeable deviations are only 
expected for energies above 300 GeV. 

In our model the residues of both soft and hard pomerons are evaluated by nonperturba- 
tive methods but this seems to agree fairly well with the phenomenological value. This has 



already been noted in ET[ for a large number of total cross sections and forward amplitudes. 



We have extended the use of our model to photoproduction of T(1S) mesons and found 
that it can yield at least a sizeable part of this supposedly hard process. 
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APPENDIX: CONDITIONS ON WAVE FUNCTIONS 

The two parameters in the wave functions of the vector mesons, N ans u, are determined 
by the normalisation condition and the leptonic decay width. 

The square of the wave functions ( |2.30| , p3l| ) summed over internal helicities is given for 
the transverse case by 

|fei(^,r)| 2 = \<Px(z,r)\ 2 x (u 4 r 2 [z 2 + (1 - z) 2 } + mfj (A.l) 
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Meson 


m f 




N 


S 




N 


S 




[GeV] 


[GeVJ 




rr i 

[fmj 


-\ t~\ 

[GeV] 




|fm| 


J/ip{lS) 


1.2 


0.59 


3.30 


0.35 


0.64 


1.47 


0.36 




1.25 


0.58 


3.17 


0.36 


0.63 


1.44 


0.35 




1.3 


0.57 


3.04 


0.36 


0.62 


1.41 


0.35 


T(1S) 


4.2 


1.29 


2.48 


0.16 


1.32 


1.17 


0.16 




4.4 


1.26 


2.35 


0.16 


1.30 


1.14 


0.16 



TABLE III: Values of u , N and size parameter S for transverse wave functions of J/tp and T 
mesons. 



and for the longitudinal case by 

|Vx,o(^r)| 2 = 2 \M^r)\ 2 x (uAzz) 2 . (A.2) 

where the X in the index stands for BSW P7| or BL -type wave functions. It has to 
fulfil the normalisation condition 

£dz J \iJxA^r)\ 2 d 2 r=l . (A.3) 

The condition which relates the wave function with the e.m. decay width fy is, for the 
transverse cases, 



M v 16tt 3 Jo 
and for the longitudinal cases 



l/ = ^^ L T^T / dz [d 2 k—^([z 2 + (l-z) 2 }k 2 + m 2 f )Mz,k) (A.4) 

My i07l s JO J Z[L — Z) v J ' 



fy = e V LuV3~^ T f 1 dz f d 2 k 16 z(l- z) 4> x (z,k) . (A.5) 

lu7T JO J 

where <px{z, k) is the Fourier transform of <fix{z, r), defined through 

r d 2 ~k ~ 

4>x(z, r) = J — 4> x (z, k) exp[-ik ■ r] , (A.6) 

and if is the quark charge in units of the elementary charge, that is e = 2/3 for the J/ip 
and -1/3 for the T. 

The decay constant fy is related to the e.m. decay width T e + e - through 



2 _ 3M v T e+e - 
tv ~ Ana" • (A - 7) 



In Table |T| we give the values of N, uo and the mean square transverse radius for the 
J/ijj and the T wave functions for different choices of the charm and bottom masses for the 
transverse BSW- and BL-type meson wave functions. 
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Meson 


m f 


C(A = ±1) 


C(A = ±1) 




[GeVJ 


[GeV J 


[GeV J 


J/ip{lS) 


1.2 


0.0103 


0.0113 




1.25 


0.0095 


0.0102 




1.3 


0.0088 


0.0093 


T(1S) 


4.2 


-0.00034 


-0.00035 




4.4 


-0.00030 


-0.00031 



TABLE IV: Values of the integrated square radius over the overlap functions, defined by Eq. 
( A.10j ), for the two kinds of transverse wave function. 



The 7*-vector meson overlap functions, necessary for the calculation of photo- and elec- 
troproduction of vector mesons are obtained from Eqs. fl2.30| , |2.3ll) and ( |2.27| , |2.28|) , and are 
given by 

a) transverse 



p y * V; i(z,r) = e v - 



'Ga 



2tt 



x (eu\[z 2 + {l-z) 2 )K 1 (er)+m 2 f K (er))=e v p jV>1 {z,r) (A.8) 



b) longitudinal 



'3a 



Prv,o(z,r) = -IGey^P 1 uj (f) X (z,r)z 2 (l - zf Q K (e r) = e v p jVi0 (z,r) , (A.9) 

Z7T 



where e = ^ zzQ 2 + m 2 

The real photon-vector meson overlap is obtained by setting Q = 0, i.e. the longitudinal 
part vanishes and e -> mj. 

The integrals of the square radius over the overlap functions 



C(X) = dz / d r r p jVjX (z,r) . 



(A.10) 



give an estimate of the effective strength for the production process, the square of these 
quantities being nearly proportional to the integrated cross sections. Thus the numbers in 
Table tell us that the ratio between T and J/ip production cross sections at 20 GeV is 
about 10" 3 . 
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